Fidelity susceptibility and long-range correlation in the Kitaev honeycomb model 
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We study exactly both the ground-state fidelity susceptibility and bond-bond correlation function 
in the Kitaev honeycomb model. Our results show that the fidelity susceptibility can be used to 
identify the topological phase transition from a gapped A phase with Abelian anyon excitations to 
a gapless B phase with non- Abelian anyon excitations. We also find that the bond-bond correlation 
function decays exponentially in the gapped phase, but algebraically in the gapless phase. For the 
former case, the correlation length is found to be l/£ = 2 sinh -1 [V2 J z — 1/(1 — Jz)], which diverges 
around the critical point J z — (l/2) + . 
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I. INTRODUCTION 

Quite recently, agreat deal of effort [H, i, S 0, S, & 0, 0, 
©, U3, El El, El El El, EH has been devoted to the role 
of fidelity, a concept borrowed from quantum information 
theory [IT]], in quantum phase transitions(QPTs)[l|j]. 
The motivation is quite obvious. Since the fidelity is a 
measure of similarity between two states, the change of 
the ground state structure around the quantum critical 
point should result in a dramatic change in the fidelity 
across the critical point. Such a fascinating prospect 
has been demonstrated in many correlated systems. For 
example, in the one-dimensional XY model, the fidelity 
shows a narrow trough at the phase transition point [2|. 
Similar properties were also found in fermionic 0] and 
bosonic systems [4|. The advantage of the fidelity is that, 
since the fidelity is a space geometrical quantity, no a 
priori knowledge of the order parameter and symmetry- 
breaking is required in studies of QPTs. 

Nevertheless, the properties of the fidelity are mainly 
determined by its leading term 0, §], i.e., its second 
derivative with respect to the driving parameter (or the 
so-called fidelity susceptibility Q). According to the 
standard perturbation method, it has been shown that 
the fidelity susceptibility actually is equivalent to the 
structure factor (fluctuation) of the driving term in the 
Hamiltonian ||. For example, if we focus on the thermal 
phase transitions and choose the temperature as the driv- 
ing parameter, the fidelity susceptibility, extracted from 
the mixed state fidelity between two thermal states [6(, is 
simply the specific heat [3, From this point of view, 
the fidelity approach to QPTs seems still to be within the 
framework of the correlation functions approach, which 
is intrinsically related to the local order parameter. 

However, some systems cannot be described in a frame- 
work built on the local order parameter. This might 
be due to the absence of preexisting symmetry in the 
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Hamiltonian, such as topological phase transitions [191 ] 
and Kosterlitz-Thouless phase transitions [2fJ. For the 
latter, since the transition is of infinite-order, it has al- 
ready been pointed out that the fidelity might fail to 
identify the phase transition point [1, EH- Therefore, it 
is an interesting issue to address the role of fidelity in 
studying the topological phase transition. 

The Kitaev honeycomb model was first introduced by 
Kitaev in search of topological order and anyonic statis- 
tics. The model is associated with a system of 1/2 spins 
which are located at the vertices of a honeycomb lattice. 
Each spin interacts with three nearest neighbor spins 
through three types of bonds, called "x(y, z)-bonds" de- 
pending on their direction. The model Hamiltonian (2lT | 
is as follows: 

H = -J x °>t~ J y E CT i °fc ~ ^ z E ^ 

aj-bonds y-bonds z-bonds 

= —JxH'x — JyHy — J Z H Z . (1) 

where j, k denote the two ends of the corresponding bond, 
and J a ,<7 a (a = x,y,z) are dimensionless coupling con- 
stants and Pauli matrices respectively. Such a model is 
rather artificial. However, its potential application in 
topological quantum computation has made it a focus of 
research in recent years 0, [H [H, [H [H, [H, H3, [M [U, 

SHU,!!!. 

The ground state of the Kitaev honeycomb model 
consists of two phases, i.e., a gapped A phase with 
Abelian anyon excitations and a gapless B phase with 
non-Abelian anyon excitations. The transition has been 
studied by various approaches. For example, it has been 
shown that a kind of long range order exists in the dual 
space (26|, such that basic concepts of Landau's theory 
of continuous phase transitions might still be applied. In 
real space, however, the spin-spin correlation functions 
vanishes rapidly with increasing distance between two 
spins. Therefore, the transition between the two phases 
is believed to be of topological type due to the absence 
of a local order parameter in real space [2l| . 

In this work, we firstly try to investigate the topolog- 
ical QPT occurring in the ground state of the Kitaev 
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honeycomb model in terms of the fidelity susceptibility. 
We find that the fidelity susceptibility can be used to 
identify the topological phase transition from a gapped 
phase with Abelian anyon excitations to gapless phase 
with non- Abelian anyon excitations. Various scaling and 
critical exponents of the fidelity susceptibility around the 
critical points are obtained through a standard finite-size 
scaling analysis. These observations from the fidelity ap- 
proach are a little surprising. Our earlier thought was 
that the fidelity susceptibility, which is a kind of struc- 
ture factor obtained by a combination of correlation func- 
tions, can hardly be related to the topological phase tran- 
sition, since the latter cannot be described by the corre- 
lation functions of local operators. So our second moti- 
vation following from the first one is to study the dom- 
inant correlation function appearing in the definition of 
the fidelity susceptibility, i.e., the bond-bond correlation 
function. We find that the correlation function decays 
algebraically in the gapless phase, but exponentially in 
the gapped phase. For the latter, the correlation length 
takes the form l/£ = 2sinh~ 1 [-\/2J 2 — 1/(1 — J z )] along 
a given evolution line. Therefore, the divergence of the 
correlation length around the critical point J z = (l/2) + 
is also a signature of the QPT. 

We organize our work as follows. In Sec. HH we in- 
troduce briefly the definition of the fidelity susceptibility 
in the Hamiltonian parameter space, then we diagonal- 
ize the Hamiltonian based on Kitaev's approaches and 
obtain the explicit forms of the Ricmann metric tensor, 
from which the fidelity susceptibility along any direction 
can be obtained. The critical and scaling behaviors of 
the fidelity susceptibility are also studied numerically. In 
Sec. IIII1 we explicitly calculate the bond-bond correla- 
tion functions in both phases. Its long range behavior 
and the correlation length in the gapped phase are stud- 
ied both analytically and numerically. Sec. HVl includes 
a brief summary. 



II. FIDELITY SUSCEPTIBILITY IN THE 
GROUND STATE 

To study the fidelity susceptibility we notice that the 
structure of the parameter space of the Hamiltonian (fT]) 
is three dimensional. In this space, we can always let the 
ground state of the Hamiltonian evolves along a certain 
path in the parameter space, i.e., 



Ja — <^a(A), 



(2) 



where A is a kind of driving parameter along the evo- 
lution line. We then extend the definition of fidelity to 
this arbitrary line in high-dimensional space. Following 
Ref. [|J , the fidelity is defined as the overlap between two 
ground states 



F = |<*o(A)|*o(A + «A)>| > 



(3) 



where <5A is the magnitude of a small displacement along 
the tangent direction at A. Then the fidelity susceptibil- 



ity along this line can be calculated as 
-21nE 



XF 



<5A^0 5X 2 



J29abn a n b , 

ab 



(4) 



where n a — dJ a /d\ denotes the tangent unit vector at 
the given point, and g a b is the Riemann metric tensor 
introduced by Zanardi, Giorda, and Cozzini 7]. For the 
present model, we have 



g a b 



E 



(tt n (A)|g a |tto(A))(tt (A)|g 6 |tt n (A)) 
{E n - E ) 2 



(5) 



where |\1/„(A)) is the eigenstate of the Hamiltonian with 
energy E n . Clearly g a b does not depend on the specific 
path along which the system evolves. However, once g a }, 
are obtained, the fidelity susceptibility is just a simple 
combination of g a b together with a unit vector which de- 
fines the direction of system evolution in the parameter 
space. 

According to Kitaev [2l|, the Hamiltonian (fT]) can be 
diagonalized exactly by introducing Majorana fermion 
operators to represent the Pauli operators as 



a x = ib x c, a y 



ib y 



c, a 



ib z c, 



(6) 



where the Majorana operators satisfy A 2 = 1, AB = 
-BA for A,B e {b x ,by,b z ,c} and A ^ B, and also 
b x b y b z c = \ip) to ensure the commutation relations 
of spin operators. Then the Hamiltonian can be written 
as 



H — UjkJ a . k Cj 



CjCfc. 



(7) 



3,k 



Since the operators Uj k — \b a - Jk b a k 3k satisfy [v,jk,H] = 0, 
[ujk,u m i] — 0, and u 2 k — 1, they can be regarded as 
generators of the Zi symmetry group. Therefore, the 
whole Hilbert space can be decomposed into common 
eigenspaces of Ujk, each subspace is characterized by a 
group of Ujk = ±1. The spin model is transformed to a 
quadratic Majorana fermionic Hamiltonian 



k CjC k . 



(8) 



Here we restrict ourselves to only the vortex free sub- 
space with translational invariants, i.e., all Ujk — 1. Af- 
ter Fourier transformation, we get the Hamiltonian of a 
unit cell in the momentum representation [21( , 



H = J2 



-1,1 

-q,2 



if (q) 
-if (q)* 



°q,i 
a q , 2 



where q = (q x ,q v ), 




(9) 



(10) 
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FIG. 1: (Color online) Fidelity susceptibility as a function 
of J z along the dashed line shown in the triangle for various 
system sizes L — 101, 303, 909. Both upper insets correspond 
to enlarged pictures of two small portions. 



r refers to the coordinate of a unit cell, 7 to a position 
type inside the cell, and 

/(q) = e q + iA q , 

e q = J x cosq x + JyCOSq y + J z , 
A q = J x smq x + JySmq y . (11) 

Here, we set L to be an odd integer, then the system size 
is N = 2L 2 . The momenta take the values 



2mr 



L-l L—l 



Qx{v) = ~' n = 2 



1 1 



(12) 



The above Hamiltonian can be rewritten using fermionic 
operators as 



H = J2 V c l + A l ( C L C ^ ~ Cla^a) ■ (13) 



Therefore, we have the ground state 

l*0> = rK. 2 io) 



= TT- , 

q V 5 ^ A q + ie 
with the ground state energy 




- q ,i + a_ q , 2 |0),(14) 



(15) 



The fidelity of the two ground states at A and A' can 
be obtained as 



A q A q + e q e q 



n^ 2 (e q -^ q ) 



FIG. 2: (Color online) Finite size scaling analysis for 
the case of power-law divergence for system sizes L = 
201, 301, . . . , 901. The fidelity susceptibility, considered as a 
function of system size and driving parameter is a function of 
L"(J Z — J™ ax ) only, and has the critical exponent u — 0.96. 



with 

cos(2fl q ) = |S-,sin(20 q ) = ^, 

-frq -& q 

The Riemann metric tensor can be expressed as 



q 

dJ b 



(17) 



(18) 



where 

6>(20 q ) 

dJ x 
d(2fl q ) 
dJ y 

d(20 q ) 
dJ z 



J z sin q x + J y sin (q x - q y ) A q 



e 2 + A 2 
q ' q 



J x sin (q x - q y ) - J z sin q y A q 



f 2 + A 2 
q ' q 



J x sin q x + J y sin q y A q 



I A c 



(19) 



Clearly, with these equations, we can in principle calcu- 
late the fidelity susceptibility along any direction in the 
parameter space according to Eq. ((4]). Here, we would 
like to point out that the same results can be obtained 
from the generalized Jordan- Wigner transformation used 
firstly by Feng, Zhang, and Xiang[26j. 

Following Kitaev [21| . we restrict our studies to the 
plane J x + J y + J z = 1 (see the large triangle in Fig. [1]). 
According to his results, the plane consists of two phases, 
i.e., a gapped A phase with Abelian anyon excitations 
and a gapless B phase with non- Abelian excitations. The 
two phases are separated by three transition lines, i.e. 
(16) Jx — 1/2; Jy — 1/2, and J z — 1/2 which form a small 
triangle in the B phase. 
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FIG. 3: (Color online) Fidelity susceptibility and a few low- 
lying excitations as a function of J z in a small portion of the 
evolution line for system size L — 51. 



FIG. 4: (Color online) Fidelity susceptibility as a function of 
J x = 2/3 — J y along the dashed line shown in the triangle for 
various system sizes L = 101, 303, 909. 



Generally, we can define an arbitrary evolution line on 
the plane. Without loss of generality, we first choose the 
line as J x = J y (see the dashed line in the triangle of Fig. 
[T]). Then the fidelity susceptibility along this line can be 
simplified as 



XF 



16 ^ 



sin q x + sin q y 



(20) 



The numerical results of different system sizes are shown 
in Fig. [T] First of all, the fidelity susceptibility per site, 
i.e. Xf/N diverges quickly with increasing system size 
around the critical point J z = 1/2. This property is sim- 
ilar to the fidelity susceptibility in other systems, such 
as the one-dimensional Ising chain [2j and the asymmet- 
ric Hubbard model [12]. Secondly, Xf/N is an intensive 
quantity in the A phase (J z > 1/2), while in the B 
phase, the fidelity susceptibility also diverges with in- 
creasing system size. Thirdly, the fidelity susceptibility 
shows many peaks in the B phase, the number of peaks 
increases linearly with the system size L (see the left 
upper inset of Fig. [1]). The phenomena of fidelity sus- 
ceptibility per site in the B phase have not been found in 
other systems previously, to our knowledge, so that they 
are rather impressive. 

To study the scaling behavior of the fidelity suscepti- 
bility around the critical point, we perform a finite-size 
scaling analysis. Since the fidelity susceptibility in the A 
phase is an intensive quantity, the fidelity susceptibility 
in the thermodynamic limit, scales as [l2j 



XF 

N 



oc 



Jz - JS 



(21) 



around J z = 1/2. Meanwhile, the maximum point of \f 
at J z = J" lax for a finite sample behaves as 



XF 

N 



(22) 



with n = 0.507± 0.0001 (see the inset of Fig. |2J. Accord- 
ing to the scaling ansatz, the rescaled fidelity suscepti- 
bility around its maximum point at J™ ax is just a simple 
function of the rescaled driving parameter, i.e., 



xT x - Xf 

XF 



f[L»(J z -J™*)}. 



(23) 



where f(x) is a universal scaling function and does not 
depend on the system size, and v is the critical expo- 
nent. The function f(x) is shown in Fig. Clearly, the 
rescaled fidelity susceptibilities of various system sizes fall 
onto a single line for a specific v = 0.96 ± 0.005. Then 
the critical exponent a can be obtained as 



a = - = 0.528 ±0.001. 
v 



(24) 



One of the most interesting observations is that a huge 
number of peaks appear in the B phase. The scaling 
analysis shows that the number of peaks is proportional 
to the system size. Physically, a peak means that the 
ground state can not adiabatically evolve from one side of 
the peak to the other side easily because the two ground 
states have distinct features. From this point of view, the 
ground state in the B phase might be stable to a adiabatic 
perturbation. Moreover, the existence of many peaks can 
also be reflected by reconstruction of the energy spectra. 
For this purpose, we choose a small portion of the evo- 
lution line and plot both the fidelity susceptibility and a 
few low-lying excitations in Fig. [3] Since the fidelity is 
inversely proportional to the energy gap [Eq. JSJ], the 
location of each peak corresponds to a gap minimum. 

Similarly we can also choose the system evolution line 
as J z — 1/3, the fidelity susceptibility then takes the 
form 



XF 



36 ^ 



(sin q x - sin q y ) + 2 sin (q x - q y ) 



A 2 
q 



(25) 
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The numerical results for this case are shown in Fig. [4] 
The results are qualitatively similar to those of previous 
cases. In the B phase, there still exist many peaks. Both 
the number and the magnitude of the peaks increase with 
the system size, while in the A phase, the fidelity suscep- 
tibility becomes an intensive quantity. 



III. LONG-RANGE CORRELATION AND 
FIDELITY SUSCEPTIBILITY 

Follow You, et al. 8], the fidelity susceptibility is a 
combination of correlation functions. Precisely, for a gen- 
eral Hamiltonian 

H = H + XH I , (26) 
the fidelity susceptibility can be calculated as 

XF = Jr [<¥ o |.Hi(T)fri(0)|¥o) - (*o|^/|*o) 2 ] dr{21) 
with t being the imaginary time and 

Therefore, the divergence of the fidelity susceptibility at 
the critical point implies the existence of a long-range 
correlation function. Without loss of generality, if we 
still restrict ourselves to the plane J x + J y + J z = 1 and 
choose J z ( J x = J y ) as the driving parameter, the bond- 
bond correlation function is defined as 

C(ri,r 2 ) = « 4 <,2<,i<, 2 ) 

-«,1<,2>«,1<,2>- (28) 

Here the subscripts ri , 1 and ri , 2 denote the two ends of 
the single z-bond at ri=(x,y). In the vortex-free case, 
through Eqs. ©, (jTOJ) , and (fT4"l) . the spin operators 
°ri i°ri 2 can be expressed in the form of fermion op- 
erators. So we finally get 



r 2 ,l a r 2 ,2 



1 N ^ E n 



(29) 



and 



(*oK,i<,2<,i<, 2 |*o) 
i^jcosKq-q'Mrr-r^-l} 

w (A q A q >-e q e q >) 

A 1 
EqEq> 



(30) 



with q q' and ri 7^ r 2 . The same results can also 
be obtained by using the Jordan- Wigner transformation 
method 0,H||. 

We show the dependence of the correlation function 
Eq. on the distance for a finite sample of L = 100 
in Fig. O Obviously, the lines can be divided into two 
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FIG. 5: (Color online) Bond-bond correlation function as a 
function of distance r for various J z and a finite sample of 
L — 100, where ri — r2 = (r, r). Downward peaks in top lines 
are due to zero-point crossing. 
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FIG. 6: (Color online) Fidelity susceptibility and the correla- 
tion function at ri — Y2 = (L/2, L/2) as a function of J z for 
a finite sample of L = 100. 



groups. In the gapless phase (J z < 1/2), the correlation 
function decays algebraically, while in the gapped phase 
(J z > 1/2), it decays exponentially. If J z < 1/2, the 
denominator in Eq. (|30|) has two zero points, which are 
of order 1/N in the large N limit. Their contribution 
causes the summation to be finite in the thermodynamic 
limit. Then using the stationary phase method, we can 
evaluate the exponents of the correlation function at long 
distance to be 4, i.e., 



C(ri,r 2 ) cx 



1 



(31) 



|rx -r 2 |* 

From Fig. the average slope of the top three lines 
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around r — 10 is estimated to be 4.05, which is slightly 
different from 4. Nevertheless, we would rather inter- 
pret the difference as due to both finite size effects and 
numerical error. On the other hand, if J z > 1/2, the 
phase is gapped and the denominator in Eq. (|30[) docs 
not have zero point on the real axis. Therefore, the whole 
summation is strongly suppressed except for the case of 
small |ri — r2 1 , whose range actually defines the correla- 
tion length. In order to evaluate the correlation, we need 
to extend the integrand (in the thermodynamic limit) in 
Eq. (|30[) to the whole complex plane, where we can find 
two singular points. Using the steepest descent method, 
we can evaluate the correlation length to be 



- = 2sinh- l V ^ . (32) 

Obviously, the correlation length becomes divergent as 
J z — ► 0.5 + . This property can also be used to signal the 
QPT occurring in the Kitaev honeycomb model in addi- 
tion to the fidelity and Chern number [2lj . The correla- 
tion length we obtained is the same as that of the string 
operators [27| . which, however, is a non-local operator. 

Although it is not easy to calculate the fidelity sus- 
ceptibility from the correlation function directly due to 
the dynamic term in Eq. (|27[) . our conjecture is con- 
firmed for the present model. That is the divergence of 
the fidelity susceptibility is related to the long-range cor- 
relations. Fig. [S] is illustrative. The correlation function 
at ri — Y2 = (£/2,L/2), in spite of its smallness, re- 
mains nonzero in the region J z < 1/2, but it vanishes in 
J z > 1/2. For the former, the oscillating structures of 
the two lines meet each other. 



sition occurrs in the honeycomb Kitaev model. Though 
no symmetry breaking exists and no local order param- 
eter in real space can be used to describe the transition, 
the fidelity susceptibility definitely can indicate the tran- 
sition point. We found that the fidelity susceptibility per 
site is an intensive quantity in the gapped phase, while 
in the gapless phase, the huge number of peaks reflects 
frequent spectral reconstruction along the evolution line. 
We also studied various scaling and critical exponents of 
the fidelity susceptibility around the critical points. 

Based on the conclusions from the fidelity, we fur- 
ther studied the bond-bond correlation function in both 
phases. We found that the bond-bond correlation func- 
tion, which plays a dominant role in the expression for the 
fidelity susceptibility, decays exponentially in the gapped 
phase, but algebraically in the gapless phase. The critical 
exponents of the correlation function in both the gapless 
and gapped phases arc calculated numerical and analyti- 
cally. Therefore, in addition to the topological properties 
of the Kitaev honeycomb model, say, the Chern num- 
ber, we found that both the fidelity susceptibility and 
the bond-bond correlation functions can be used to wit- 
ness the QPT in the model. 

Note added. After finishing this work, we noticed that 
a work on the fidelity per site instead of the fidelity sus- 
ceptibility in a similar model appeared [33|. 
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